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We explore the consequences of Replica Symmetry Breaking at zero temperature. We 
introduce a repulsive coupling between a system and its unperturbed ground state. In the 
Replica Symmetry Breaking scenario a finite coupling induces a non trivial overlap probabil- 
ity distribution among the unperturbed ground state and the one in presence of the coupling. 
We find a closed formula for this probability for arbitrary ultrametric trees, in terms of the 
parameters defining the tree. The same probability is computed in numerical simulations of 
a simple model with many ground states, but no ultrametricity: polymers in random media 
in 1+1 dimension. This gives us an idea of what violation of our formula can be expected in 
cases when ultrametricity does not hold. 



In recent times there has been a wide interest in the behavior spin glasses with Gaussian couplings at zero 
temperature. 

Some of the reasons for this interest are the following: 

• The energies are continuous variables and the ground state is unique. It is also natural to suppose 
(although it is far from being proved) that the limits T — > and N — > oo do commute and therefore the 
shape of the energy landscape is similar to that of the free energy landscape at non-zero temperature 
(for a discussion of this point see |Q). 

• Working at zero temperature avoids completely the possibility that the temperature used is too near 
to the critical point. 

• Technical progresses has been done in the algorithm for finding the ground state |^|J||||7]|| and it is 
now possible to studies three dimensional systems up to 14 3 spins Q. 

In this framework it has been suggested that a possible test of the applicability of the Replica Symmetry 
Breaking (RSB) scenario is the study of the overlap of the ground state of two systems whose total Hamiltonian 
differs by a quantity of order 1 MM . 

Let us consider a simple case. We have a first system with Hamiltonian Hq (er) and its ground state is given 
by Tj. We now consider a second system whose Hamiltonian is 



I. INTRODUCTION 



H x {a) = Ho(a) + eH T (a) . 



(1) 



Three quite simple choices of H T (a) are: 



H r {o) = q{cr,T), H T (cr) = q 2 (a,T 



), H T (a) = qi{cr,r) , 



(2) 



where the overlap q and the link overlap qi are given by 





(3) 



(i,k) 
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where the sum is done over all the nearest neighbor pairs (i, k) in a short range model or over all the pairs in 
the SK model (Ni being the total number of pairs in this sum) . The third possibility has been actually used 

by|. 

In presence of quenched disorder, the the value of the overlap among the ground states of Hq and H\ can 
be sample dependent. This observation can be used tp the starting point for investigating possible RSB in 
the three dimensional Edwards- Anderson model. The question we address in this paper, is the computation 
of the probability distribution induced by the random couplings of q or of qi among the two ground states, 
in the hierarchical RSB scenario. 

Obviously the choice H q (a) — q is interesting only in presence of a magnetic field which breaks the 
symmetry a(i) — > —a(i), otherwise we would get that a(i) = —t(i) for positive non zero e and q = — 1. The 
second choice is more interesting at zero magnetic field, but it is slightly harder to implement numerically, 
because its non local nature. The third choice is however equivalent to the first one in the SK model, where 
is known that qi = q 2 apart from corrections that vanishes when the number of spins goes to infinity. In short 
range models, it is possible (as suggested by the principle of replica equivalence [§|) that with probability one 
when the volume goes to infinity qi = f{q 2 ), where / is a function that can be determined numerically and 
which should be not too far from 



f(q 2 ) = A+(l-A)q 2 . (4) 

It is evident that for finite e the perturbation is of order 1 and it quite interesting that if replica symmetry is 
broken the function P(q) is non trivial at e ^ 0. Let us define as E gs (0) + A(g) the energy of the first excited 
state of the Hamiltonian Hq with an overlap q with the ground state. The ground state of the Hamiltonian 
H x is 

E gs (e) = E gs (0) + min {A(q) + eg(q)} (5) 

0<<j<1 

where g(q) = q,q 2 ,qi. The main achievement of this paper will be the computation in section II and III of 
the joint probability distribution of A and q for which the minimum is attained, for arbitrary RSB trees. 
In section IV we will give some example in mean field models, while in section V we show the result of a 
numerical computation for directed polymers in random media in 1+1 dimension. 



II. REPLICA SYMMETRY BREAKING 



The computation presented in this note could be done in two different ways: 

• Using the replica formalism |^] to compute the partition function of the perturbed Hamiltonian Hi . 

• Exploiting directly the information coming from replica symmetry breaking on the probability distri- 
bution of the lowest lying states and doing a pure probabilistic computation. 

The first alternative leads to some apparently messy combinatorial analysis so that we have decided to follow 
the second alternative. In this case the computation is physically instructive. 

In this section we will recall, without proof, some known results about replica symmetry breaking and also 
find some new consequences of those results. Let us assume that replica symmetry is broken in the system 
we consider and its breaking is characterized by a function x(q,T) such that in the low temperature limit 

x(q,T) = Ty(q) + 0(T 2 ) , (6) 

where the function y(q) may be singular at q = 1 (in the SK model it diverges as (1 — q)^ 1 ^ 2 near q — 1) ]l0| . 
The space of lower lying configurations is organized in a rather complex way. 



A. One step replica symmetry breaking 

In this case only two values of the overlap are allowed (qo and 1), i.e. all different minima have a mutual 
overlap equal to qo- If we call R a reference total energy, which depends on the choice of the systems, i.e. on 
the variables J, the probability to find a configuration in the interval (E, E + dE) is given by 
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vq{E\R) = exp(yo(E - R)) . (7) 

We notice that configurations which differs by a number of spin flips which remains finite when the volume 
goes to infinity are identified. 

This well known results has the consequence that the probability distribution of the ground state energy 
Eq is given by the Gumbel law 

Ho(E \R) = exp(y (£o - R)) exp(-A exp(y (E - R)) , (8) 

with A = 1/j/o- This formula is easily understood noticing that the probability that there are no configura- 
tions for E' < E is given by 

exp ^- J dE' eMvo(E' - R)^j = cxp(-y- 1 exp{y Q (E - R)) . (9) 

In the same way we obtaining that the probability of having a ground state at E and the first excited 
configuration at E\ is given by: 

P{Eq,E x \R) = is (E \R)MEi\R) , 

= cxp(y (£o - R)) eMvo(Ei - R)) exp^ 1 exp(2/o(^i - R)) ■ (10) 

Finally if we define A = E\ — Eq, the probability distribution of A , integrated over E and E\, one finds 
the simple result 

P(A ) =y - 1 exp(-y Ao) , (11) 
for positive A , the probability being obviously zero for negative Ao. 

B. Many level replica symmetry breaking 

Let us in this section generalize the computation for an arbitrary RSB tree (see figure). As customary, we 
will first consider a tree with k levels and at the end we will generalize the result to the continuous branching 
limit. The construction of the tree has been described many times, and we only repeat it briefly to fix the 
notation. At each node of the tree at the level I it is assigned an energy Ei, which is chosen in such a way 
that the number of nodes with energy in the interval (Ei, Ei + dE{) branching from a node with energy Ei_i 
is a Poisson variable with average equal to exp(yi(Ei — Ei_\))dEi. We will consider the case in which for all I, 
Vi+i > Vi, which will be a necessary condition of convergence of the integrals that appear in the computation. 
The root energy E = R is the reference energy of the previous section. We call ^-clusters, the set of branches 
which coincide at the Z-th level of the tree. 




E_gs E_gs+ A_2 E_gs+A_l E_gs+A_0 

FIG. 1. A three level tree, expliciting the notation we use in the text. 



Let us define for any I = k — 1, 1 the first ^-excited state, as the first excited state which is in the same 
Z-cluster as the ground state, but in a different I + 1-cluster, and denote its energy E gs + A;. In this section we 
compute the joint probability distribution of all the ^-gaps A/. We will get this quantity by first computing 
Pk(E gs , E gs + Afc, E gs + Ai\Eq) and then integrating over E gs . In that computation we make use of the 
following properties: 

1. The probability of a ground state fj,i(E\Ei) of an i-cluster (we call it an I ground state) is given by 

m{E\E t ) = exp (y l+1 (E - E[) + A ie y^ E - E 'A . (12) 



where the A\ are a positive constants whose value could be easily compute, but we will not need. For 
I = k the formula was derived in the previous section. Let us now proceed by induction supposing that 
the formula holds for I + 1 and show that it holds for I. Under the induction hypothesis, we find that 
the number of I + 1-ground states with energy E in an I cluster is given by 

vi +1 (E\Ei)dEi = J dE l+1 e m +^ El +i- E ^ IM+1 (E\E l+1 ) 

= const x e yi + l( - E - El) dEi (13) 



from which we immediately find that the distribution of the I ground state is given by (|12j) exploiting 
the reasonings of the previous sub-section. 

2. The joint probability Pk{Eg Sl E g s + A k , ...,E gs + Ai\E k , Ek-i, E ) can be written as: 

Pk(E gs , Egs + Afc_i, ...,Eg S + Aol-Bfe-i, ...,-Eo) = 

i/ fe (£ ss |£fc_iK(£ gs + A fe _ 1 |^ fe _i) M/s _i(^ gs + A fe _ 2 |£fc-2) x ... x m(E ga + A \E Q ) (14) 
from which we get: 
Pk{E gs ,Egs + Afc, Egs + AoI-Eq) = 

J &E k _ l ..AE 1 v k (E ga \E k - 1 )n k (E ga + A k \E k _ 1 )e Vk -^ E "-'- Ek -^ l i k - 1 {E gs + A k _ 2 \E k _ 2 )e^-^ E ^- E "-^ x 
... x u 2 (E gs + A 2 \E 1 )ey^- E ^n 1 (E gs + Ai|£„) (15) 

3. A detailed computation shows that 

JdE^ u k (E gs \E k ^)tx k (E gs + A k \E k _ 1 )e y *-^ El '-i- E >°-^ = (y fc - y fc _ 1 )e-^-^- 1 ) A ^ fc _ 1 (^|-B fe _ 2 ). 

(16) 

This allows to integrate all the Ei (I = k — 1, 1) telescopically, and obtain 

k 

P k {Eg S ,Eg S + A fe , ...,E gs + A j_ \ E ) = n (( W - y^e-^-y-^ Vl (E gs \E )m{E ga + A^Eq) (17) 

z=2 

4. We can finally integrate E gs as in the previous section and get that the gaps' probability distribution 
is: 

k 

P(A k , Ai) = [] ((Vi - yi-i)e-^-y^ Ai ) (18) 
t=i 

having defined yo = 0. 

We can now consider the continuum branching limit, in which the branches can be indexed by the value 
of q, or by any monotonically increasing function of q. yi — > y(q), yi — Ui-i —* 2/'(<z) dq where the previous 
formula reduces to 



P({A(q)}) = cxp (- £ dqy / (q)A(q^j f[ y> \q) dq , 



(19) 

where, if y\ — > y(0) ^ 0, we make the convention that y'(0) = y(0)S(q). 
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III. THE OVERLAP PROBABILITY DISTRIBUTION 



We are finally in the position to compute the joint distribution of the gap and of the overlap. For a given 
sample, as we said, the difference among the ground state energy of the Hamiltonians Hi and Hq is given by 
the 

E gs {e) - E gs (0) = min A(q) + eg(q). (20) 

0<g<l 

Noticing that the indexing of the tree could be done by the function g(q) itself, we can concentrate here to 
case g(q) = |g|, and consider only positive overlaps. All the other cases can be obtained by this one via a 
simple change of variable. Let us call A and q the arguments of the minimum in (|2C|). Notice that if we want 
q 7^ 1, A has to verify the inequality A < e(l — q), which expresses the fact that A(l) is always equal to zero. 

We get the probability P(A,q) integrating in formula ( |l9|) all the A(q') (q' ^ q) with the condition 
A(q') > min{0, A + e(q — q')}, so that, if q' > A/e + q the integration over A(q') will contribute with a one, 
while, in the opposite case q' < A/e + q it contributes with the factor exp(— edq' y'{q'){A/e + q — q')). We 
finally get: 



P(A, g) = 0(1 - q - A/e)y'(q) exp (-e jf 



dq' (y(q') - y(0)) + d(A)S(q ~ 1) exp(-e X ) , (21) 



where we defined \ = y(0) + J Q dq y(q). The factor y'(q) that multiply the exponential comes from the only 
A which has remained unintegrated. 

Let us notice that the formula depends on A only in the combination A/e + q. 

Integrating over A we get, for the overlap probability 



P(q) = S(q - 1) exp(- ex ) + ey'(q) / dq' exp -e / dq"(y(q") - y(0)) (22) 
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It is also interesting to study the probability distribution of w = A/e + q. Integrating over A and q for 
fixed w with the condition < q < w we find the remarkable formula: 

P(w) =6(1- w)e(y(w) - y(0)) exp (-e J dq (y(q) - y(0)fj + 6(w - 1) exp(-e X ) . (23) 

The primitive of this function has a very simple dependence on y and e. If we define Q(w) — J w dw' P(w') 
we find: 

Q(w) = exp (-e £ dq (y(q) - y{0))j (24) 

which is particularly well suited for the extraction of the function y(q) from numerical simulations. Notice 
that eq. d24| ) implies that e _1 ln(Q(u')) is e-independent, which is a non-trivial result. 

IV. EXAMPLES 

In this section we show how our formula looks like in some cases. The first example we make is the one 
of spherical models. These models are defined by a random Gaussian Hamiltonian H(a) (a = {a±, ...,<7at}), 
which have correlation function 



H(a)H(r)=Nf(q(a,T)) (25) 

and the spins are subject to the spherical constraint J2i a 1 — N . In these models, the function y(q) is 
temperature independent and equal to |]llf 

y(q) = l/2f'"( q )/(f"(q)^) . (26) 



5 



This last equation makes sense only if the resulting function y{q) is an increasing function of q. This in 
particular happens for f(q) = l/2(q 2 + aq p ) if p > 4 and a small enough, where we find 



y(q) = 



a 1 



r,P-3 



3/2 



(27) 



while 



Y(w) = / dqy(q) = 1 



1 + 



(iqP 



(28) 



2(p-2)\ 



In figure |^ we show the function Y(w) and the function Q(w) for various values of e in the case p — 5, 
a = 0.3. In figure || we show the function P(q) for the same values of the parameters. 




FIG. 2. The function Y(w) in the spherical model with p — 5 a — 0.3. In the inset, the function Q(w) for e = 2 k , 
from top to bottom k = 4, 7, 10, 13, 16 



0.2 0.3 



0.6 0.7 0.8 0.9 



FIG. 3. The function P(q) for the same model and parameters of figure ^. 

In the case of the SK model the function y(q) at low temperature, has be estimated in |l0) using the 
so-called PaT approximation, and displays a square root divergence at q = 1, while starting linearly at q = 
(a best fit of the form y(q) = gives a = 1.309, b = -0.695). 

Using that estimate we immediately compute the function Q(w) which is plotted in figure ^ for various 
values of e. 



G 



FIG. 4. The function Y(w) for the SK model in the PaT approximation. In the inset, the function Q(w) computed 
by Q for e = 0.01 x 2 k with (from top to bottom) k = 3, 5, 7, 9, 11, 13. 



It is interesting to study the limit e — > oo of our formula. Let us suppose that y behaves as y(q) — aq a + .... 
for a > 1 for low q. In this case P(q) will be dominates by the behavior of y close to q = 0. We can introduce 
a cut-off A such that A — > and A 1+Q e — > oo. 



,A 

P(g) ~ eag"" 1 / exp 

Jq 



n.a+1 



a + 1 

which, rescaling the integration variable and sending the cut-off to zero becomes: 



(29) 



1 f eaq a+1 \ 1+ " 1 



at 



q \ 1 + a J 1 + a a + 1 

where the incomplete gamma function is 

poo 

T(n,x)= / dyy^e-y. (31) 

J X 

This case is relevant in the spherical models where y(q) ~ 4( p -3)! 9 P ~ 3 an d f° r the SK model where y(q) ~ ag. 

V. DIRECTED POLYMERS IN 1 + 1 DIMENSION 

The natural play ground of the exposed theory are finite dimensional s pin gl asses. Our analysis predicts 
a very peculiar dependence of the probability of q and w on e. Formulae (|22| , |24| ) has been derived assuming 
that the low l ying s tates verify ultrametricity. It is therefore interesting understand what violations of the 



scaling forms (22 2J) can be expected when the ground states structure is nontrivial, but ultrametricity does 
not hold. 

In this paper we study numerically the simple case of directed polymers in random media in 1+1 dimension. 
The model we use is defined on the square lattice, where the polymer can perform a random walk starting 
from the origin. On each site of the lattice is defined an passage energy cost which is a Gaussian variable 
with unit variance, and independent of all the other energies. The properties of this sort of models have been 
studied extensively |fl2|| , and it is well known that while the low temperature thermodynamics is dominated 
by a single ground state, there exist many "pure states" (i.e. metastable states separated by growing barriers) 
with typical energy gap with the ground state scaling as L 1 / 3 (for a polymer of total length L). The overlap 
for two polymers of length L with a common source in the origin is often defined as the fraction of the 
monomers passing in the same sites of the lattice in the two polymers. Given the previously mentioned 
scaling of the energy gap it is natural that the scaling of the coupling in order to have a non trivial P(q) as 
defined in the previous section is 

e = V L 1/3 . (32) 



If we had to suppose the validity of the formula (|22j) we would conclude that the function y(q) for samples 
of length L depends on L and scales as L -1 / 3 . 
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As we stressed, ultrametricity does not hold in this model. Although we did not make a systematic study 
we can easily show the lack of ultrametricity generalizing the coupling procedure to a third "replica" , which 
has a repulsion both with the unperturbed ground state, and with the one obtained with the coupling. For 
simplicity, we look at the case in which all couplings are equal. In figure || we show the overlap between the 
second and the third replica as a function of the overlap between the first and the third, fixing the overlap 
between the first and the second to 0.8, and we see no sign of ultrametricity. 




0.2 0.4 0.6 0.8 1 

FIG. 5. We plot here 523 as a function of 513 for q\2 = 0.8 for 1151 different samples with L — 200 and e = 0.01. It 
is apparent that ultrametricity is violated; if it was obeyed the points would stick on the line ^23 = 913 for qis < 0.8, 
on the line 923 = 0.8 for §13 > 0.8, and on the line qi3 = 0.8 for 523 > 0.8. The lines 523 = qia ± 0.2 represent 
the bounds imposed by the triangular inequality. Data obtained for L = 400 show that there is no trend towards 
ultrametricity as L is increased. 

In figure ^ we show the function Q(w) for various values of L and e of the form (^2|). We see that for values 
of L ~ 200 the expected independence of Q(w) of L is reasonably obeyed. A close inspection to Q(l) however, 
which represents the probability of q = 1 reveals that this quantity behaves as Q(l) = exp(— r]xL 2 ^ 3 ), with 
X = 0.85 ± 0.02, and that the scaling is violated in proximity of w = 1. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.6 0.9 1 



FIG. 6. The cumulative probability Q(w) and 5 different values of 77 = 0.1 x 1.4™, from top to bottom 
n = 1,5,10,15,20 and L = 100 (dotted line) L = 200 (dashed line) L = 400 (full line). The data are obtained 
on sets of 10000 different samples. We see that the scaling with L is reasonably obeyed. 

We next try to scale the data with e according to the formula (|24|) . In figure we see that the works quite 
well for small values of e, while we show in figure H that there are important violations for large values of e. 
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10000 




0.1 0.2 0.3 0.4 0.5 0.6 07 0.8 0.9 1 

FIG. 7. The cumulative probability —l/eLog(Q(w)) and 7 different values of r\ = 0.1 x 1.4™, n = 1, 7 and 
L — 400. The data are obtained on sets of 10000 different samples. We see that for these low values of e the scaling 
is reasonably obeyed. In the inset we plot Q(w) for the same values of the parameters. 




0.001 1 1 — 1 — 1 — 1 — ' ' — ' ' — 

0.0001 0.001 0.01 0.1 1 

FIG. 8. The cumulative probability — e _1 ln(Q(w)) and 5 different values of r\ = 0.1 x 1.4™, n = 15, ...20 and 
L — 400. The data are obtained on sets of 10000 different samples. In the inset the function Q(w). For these values 
of e the scaling form of the function Q(w) is violated. 



VI. CONCLUSIONS 

Monte-Carlo simulations of finite dimensional spin glasses, show a behavior in agreement with RSB fill . 
However, the use of Monte-Carlo has been criticized on the ground that one can only equilibrate the system 
too close to the critical point, where finite size effects are large and could spoil the conclusions about ergodicity 
breaking in the thermodynamic limit Jlifl (see however Jl^Jl^j ). It is important therefore to find consequences 
of RSB at zero temperature. In this paper we have devised some of them. 

We have found that a universal formula holds for the probability of the overlap among the uncoupled 
ground state and the coupled one. The investigation of the validity of that formula in three dimensional 
systems is not beyond reach with the present technology, and will furnish an important test to understanding 
the nature of the spin-glass phase of three dimensional systems. 

Two main ingredient will be involved in the calculation: the exponential distribution of the states Jl7t 
and their ultrametric organization |l^] . This implies that the tree of states is described by a single function 
y(q). The function y(q) may in principle depend on N. Mean field theory predicts that y(q) remains 
finite in the thermodynamic limit implying that the energy differences between pure states remain finite in 
the thermodynamic limit. However, one could envisage systems where both exponential distribution and 
ultrametricity are valid, but the typical energy differences scale as L . In this case, one still have a function 
y(q) which scales as L~ e and in order to measure a nontrivial overlap distribution one needs a coupling of 
order L . The numerical study of the overlap and gap distribution, and the comparison with the formulae 
found in this paper will give important information about the organization of the states in short range models. 
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